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Abstract. X-ray observations are a direct diagnostic of fast electrons produced in solar flares, 
energized during the energy release process and directed towards the Sun. Since the properties 
of accelerated electrons can be substantially changed during their transport and interaction 
with the background plasma, a model must ultimately be applied to X-ray observations in 
order to understand the mechanism responsible for their acceleration. A cold thick target 
model is ubiquitously used for this task, since it provides a simple analytic relationship between 
the accelerated electron spectrum and the emitting electron spectrum in the X-ray source, 
with the latter quantity readily obtained from X-ray observations. However, such a model 
is inappropriate for the majority of solar flares in which the electrons propagate in a hot 
megaKelvin plasma, because it does not take into account the physics of thermalization of 
fast electrons. The use of a more realistic model, properly accounting for the properties of the 
background plasma, and the collisional diffusion and thermalization of electrons, can alleviate 
or even remove many of the traditional problems associated with the cold thick target model 
and the deduction of the accelerated electron spectrum from X-ray spectroscopy, such as the 
number problem and the need to impose an ad hoc low energy cut-off. 


1. Introduction and the observation of solar flare X-rays 

During a solar flare, a very large number of electrons, of the order of 10^® electrons s“^ are 
accelerated, but the mechanism responsible for their acceleration remains poorly understood 
(see [T] for a recent review). The properties of accelerated electrons are mainly deduced by X- 
ray observations, currently using X-ray imaging and spectroscopy provided by the Ramaty High 
Energy Solar Spectroscopic Imager (RHESSI) [2]. Solar flare X-rays are mainly produced as 
bremsstrahlung from electron-ion collisions and are optically thin in the solar atmosphere where 
hydrogen number densities range from around 10® — 10^® cm~®. X-ray observations thus provide 
an untainted inference of the emitting electron distribution, which is generally different from the 
accelerated electron distribution, since the latter can be modified en route by many processes in 
the solar atmosphere. These processes include Coulomb collisions (mainly with other electrons), 
and also non-collisional and turbulent processes (for example OH]). Apart from rare cases, such 
as the ‘cold flare’ [5], most solar flare spectra consist of two parts, usually treated as separate 
entities: a bremsstrahlung power law component dominant above ~25 keV that is produced by 
the collision and deceleration of fast accelerated electrons, and a thermal component (mostly 
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Figure 1. Left: RHESSI X-ray images of two solar flares, where the top image shows a 
‘standard’ flare type (SOL2013-05-13T16:04) with a hot coronal source observed at energies of 
~ 10 to 25 keV and X-ray footpoints above ~ 30 keV located lower in the atmosphere. The 
bottom image shows a coronal dominated flare type (SOL2005-08-23T14:30), where both the 
thermal and nonthermal emission exist cospatially. Right: RHESSI X-ray spectra for each flare. 
Most flares, independent of the flare type have two main components: a hot thermal component 
dominant at lower energies (red) and a power law nonthermal component at higher energies 
(blue). 


consisting of bremsstrahlung, but also free-bound and line emission), produced by a highly 
heated plasma (1 keV or more) that completely dominates the X-ray spectrum at energies 
below 25 keV (see Figured]) [6l[7|. RHESSI X-ray imaging usually shows that the bulk of 
hot material ~ 20 MK (~ 1.7 keV) emanates from a coronal source and the nonthermal X-rays 
from ‘footpoint’ sources, located lower in the atmosphere (see Figure [T|)- There is also some 
evidence of hot emission of around 10 MK (~ 0.86 keV) closer to the footpoints [8]. Sometimes 
the thermal and nonthermal components are cospatial and are both produced in the hot, dense 
coronal loop, e.g. Pdoidl] (see Figure dj)- 

































Figure 2. Equations ([I|) and m plotted for left: a power law {nVF){E), and right: a 
Maxwellian plus power law {nVF)(E), where blue: frictional first order term (Equation (j2]) ). 
red: diffusive second order term (Equation ([2|)) and black: both terms (Equation ([2])). The 
green curve represents the cold model (Equation ([I|)). In the full model, the diffusive term at 
low energies cancels the frictional term providing a natural cutoff to the accelerated electron 
distribution, removing the low energy cutoff problem. 


2. Finding the flare accelerated electron distribntion 

Commonly, either by a forward fitting or inversion technique, the cold thick target model e.g. 
[la [la E] is applied to RHESSI X-ray spectra. This model is popular due to its simple 
analytical form. It assumes that: 1. the accelerated electrons are transported through a ‘cold’ 
material, whereby the accelerated electron energy E ^ T, where T is the temperature of the 
background material in keV, and 2. the observed X-ray bremsstrahlung emission is produced 
in a region of sufficiently high density that the accelerated electrons lose all of their energy 
collisionally by friction in this region. In this model, the dynamics of X-ray emitting electrons 
is dominated by deterministic collisional energy loss and, hence, the spatially averaged electron 
spectrum {nVF) [15] is related to the accelerated electron spectrum Eq injected into the source 
via. 


Fo{Eo) 


K d UnVF){Ey 
aIe [ E . 
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where E and Eq are the observed and injected electron energies respectively, K is the collisional 
parameter and A is the injection area. Equation ([T|) is plotted in Figure [2] for two different 
{nVE). Solar flare X-ray spectra, in the nonthermal domain, are typically quite steep (Figure 
[T|). Using the above model that only retains the effect of deterministic energy loss requires that 
the injected electron flux spectra Eq{Eq) are similarly steep. Since the total injected power is a 
diverging quantity for such steep power-laws, the concept of a ‘low-energy cutoff’ is frequently 
used to keep the total number and power of the accelerated electron distribution hnite. Such 
a quantity is also difficult to determine from observation due to the large thermal component 
dominant at lower energies (Figured]). 


2.1. Hot plasma in the corona 

The cold thick target model is not directly applicable to the analysis of hot flaring plasma, since 
it does not account for the different behaviour and thermalization of electrons at energies close to 
F ~ T. Such electrons are dominated by collisional diffusion and not by systematic energy loss, 
thermalizing with the background plasma. In order to improve the estimation of the accelerated 
electron properties, we must take into account the presence of a finite target temperature and 































Figure 3. The resulting spatial distributions (parallel to a guiding magnetic field) and spatially 
integrated spectra for a simulated accelerated electron distribution (top panel) transported 
through two different solar flare coronal environments (middle panels) of left: T = 20 MK and 
n = 10^® cm“^, and right: T = 20 MK and n = 10^^ cm“^, and injected across 10'^ at the loop 
top apex {z = 0''). The grey regions denote a low temperature, high density ‘chromospheric’ 
region. High coronal densities lead to the bulk of electrons stopping in the corona (coronal 
thick-target flares), but a low density leads to the bulk of electrons stopping in the dense low 
atmosphere (footpoint dominated flares). Bottom panels: the resulting X-ray spectrum (pink) 
is compared with that of the cold thick target model (blue). The form of the spatial distribution 
and the integrated X-ray spectrum is dependent on the properties of the background plasma 
and cannot be ignored if a more correct form for the accelerated electron distribution is to be 
determined from observation. 


























the resultant collisional diffusion of electrons. In such a model, the observed electron spectrum 
{nVF){E) and the accelerated electron spectrum Fq are related by 


Fo{Eo) — 
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where u = y/E/T, G is known as the Chandrasekhar function (see [16] for more details) and 
temperature T is in units of keV. Equation ([2|) reduces to the cold target form, in the high 
energy, low temperature limit when E ^ T, and has two very important distinctions compared 
to the cold target equation (Equation ([U). Eirst, the first order friction term is modified so that 
electrons can either lose or gain energy during a collision, which is particularly important at 
E <T (see m), and even more importantly, there is the inclusion of a second order term that 
describes energy changes by collisional diffusion, controlling the behaviour of electrons at E ~ T, 
and ultimately leading to thermalization. Equations m and ([2]) are shown in Eigure[2]for two 
different {nVF){E): a power law and a Maxwellian plus power law. Importantly, Eigure[2]shows 
that the diffusive term cancels the deterministic term at low E. Eor the case of a power law 
{nVF){E) extending over all E, the injected electron flux Fq{Eq) becomes unphysically negative 
at low injected energies Eq, producing a natural cutoff to the injected distribution. In the more 
plausible ‘Maxwellian plus power law’ case (as might be expected in a hot target). Equation (l2|) 
disappears for the Maxwellian component, once again removing the need to apply an arbitrary 
low energy cutoff to the accelerated electron distribution. 


2.2. Changes to the observed electron spectrum due to the presence of hot flaring plasma 
In order to show how the presence of hot plasma in a flaring solar atmosphere changes the 
form of the accelerated electron distribution, Eigure [3] plots the resulting spatial distributions 
and spatially averaged spectra for the same injected electron distribution, for the same target 
temperature T = 1.7 keV (20 MK) but with number densities of a) n = 10^*^ cm“^ and b) 
n = 10^^ cm“^. Unlike in the cold thick target model, the thermalization of fast electrons 
produces a thermal component at energies <10 keV that increases with the background number 
density n. Increasing the background temperature, for a given density n, will cause the thermal 
component to be more pronounced at higher energies (see |16]1. Eigure [3] exemplifies how the 
background properties change the ‘collisionaly relaxed’ electron distribution in the target. Erom 
Eigure O it is easy to recognize the thermalized part in {nVF)[E) at low energies. Importantly, 
we find that the more realistic model with diffusion correctly describes the electron energy loss 
rate, particularly at E ~ T. Hence, fewer electrons are required to produce a given, observed 
X-ray flux, of an order of a magnitude (see [ISj). The number of electrons required to produce 
a given X-ray flux falls with increasing temperature T and number density n. 


3. Discussion and ongoing work 

The work summarized demonstrates that a more realistic model that includes the effects of 
collisional diffusion and thermalization of fast electrons, is preferable for inferring the properties 
of solar flare accelerated electrons at low energies of 10 to 30 keV from X-ray observations. 
Accounting for such effects can remove problems such as the low energy cutoff, and reduce 
the number of electrons required to produce a given observed X-ray flux. Importantly, such 
a model can be applied to all flare types showing the presence of hot plasma, whether they 
are ‘standard’ footpoint-dominated or coronal-source-dominated flares. A full description of the 
work is provided in |18] . and also in [16| . Eurther, the application of such a model to real X-ray 
flare data is currently the topic of ongoing work. We hope it will eventually replace the use of the 
cold thick target model to deduce the accelerated electron distribution from X-ray spectroscopy. 






3.1. Acknowledgments 

NLSJ was funded by STFC and a SUPA scholarship. EPK and NHB gratefully acknowledge 
the financial support by an STFC Grant. AGE was supported by NASA Grant NNX14AK56G. 
NLSJ wishes to thank the organizers for providing her with the opportunity to discuss this work 
at the conference. 

4. References 

[1] Zharkova V V, Arzner K, Benz A O, Browning P, Dauphin C, Emslie A G, Fletcher L, Kontar E P, Mann G, 

Onofri M, Petrosian V, Turkmani R, Vilmer N and Vlahos L 2011 Space Sci. Rev. 159 357-420 (Preprint 
1110.2359) 

[2] Lin R P, Dennis B R, Hurford G J et al. 2002 Solar Phys. 210 3-32 

[3] Kontar E P, Bian N H, Emslie A G and Vilmer N 2014 Astrophys. J. 780 176 (Preprint 1312.0266) 

[4] Kontar E P, Ratcliffe H and Bian N H 2012 Astron. Astrophys. 539 A43 (Preprint 1112.4448) 

[5] Fleishman G D, Kontar E P, Nita G M and Gary D E 2011 Astrophys. J. Lett. 731 L19 (Preprint 1103.2705) 

[6] Kontar E P, Brown J C, Emslie A G, Hajdas W, Holman G D, Hurford G J, Kasparova J, Mallik P G V, 

Massone A M, McGonnell M L, Piana M, Prato M, Schmahl E J and Suarez-Garcia E 2011 Space Sci. 
Rev. 159 301-355 (Preprint 1110.1755) 

[7] Holman G D, Aschwanden M J, Aurass H, Battaglia M, Grigis P G, Kontar E P, Liu W, Saint-Hilaire P and 

Zharkova V V 2011 Space Sci. Rev. 159 107-166 (Preprint 1109.6496) 

[8] Graham D R, Hannah I G, Fletcher L and Milligan R O 2013 Astrophys. J. 767 83 (Preprint 1302.2514) 

[9] Xu Y, Emslie A G and Hurford G J 2008 Astrophys. J. 673 576-585 

[10] Kontar E P, Hannah I G and Bian N H 2011 Astrophys. J. Lett. 730 L22 (Preprint 1102.3664) 

[11] Guo J, Emslie A G and Piana M 2013 Astrophys. J. 766 28 (Preprint 1303.1077) 

[12] Brown J G 1971 Solar Phys. 18 489-502 

[13] Syrovatskii S I and Shmeleva O P 1972 Sov. Astronom. 16 273 

[14] Emslie A G 1978 Astrophys. J. 224 241-246 

[15] Brown J G, Emslie A G and Kontar E P 2003 Astrophys. J. Lett. 595 L115-L117 

[16] Jeffrey N L S, Kontar E P, Bian N H and Emslie A G 2014 Astrophys. J. 787 86 (Preprint 1404.1962) 

[17] Emslie A G 2003 Astrophys. J. Lett. 595 L119-L121 

[18] Kontar E P, Jeffrey N L S, Emslie A G and Bian N H 2015 ArXiv e-prints (Preprint 1505.03733) 


